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Abstract 

Using Jacobi's identity, a simple formula expressing Bessel functions of integer order as simple 
combinations of powers and hyperbolic functions, plus higher order corrections, is obtained. 
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In this article we shall propose a simple formula expressing the modified Bessel functions 
of integer order, I n , in terms of powers and hyperbolic functions of the same argument. It 
can be easily adapted for the Bessel functions J n . 

The starting point is the generalization of the Jacobi identity ([l| p. 22) used in the 
calculation of lattice sums [2| : 
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For w — 1 and N — 2, (JI1) gives a well-known formula: 
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COsh Z = ^i Z ) = W + 2 [ J 2(^) + + , 

fc=— 00 

(see for instance [3(, eq.9.6.39) 

It is interesting to exploit (jTJ) at w = 1, for larger values of N. For = 4, 
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and for N = 8 
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(4) 



-^)=I (s) + 2 

It is easy to see that, if the l.h.s. of these equations contains p hyperbolic cosines, it 
provides an exact expression for the series of Iq(z), cut off at the z 4p term. The generalization 
of (jlj) for TV = 4p, with p - an arbitrary integer, is indeed: 
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Because 



IJ-iz) = i n JJz). 



I (z) + 2j2h P k(z). 
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our result (ED can be written as: 
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It is easy to obtain formulae similar to (jSJ) for any modified Bessel function of integer 
order. Let us introduce the notations: 

71 p — 1 

Ci=COS— , ... C p _i = COS— 7T. (8) 

2p 2p 

and let us define the functions: 



S q (z) = sinh z + 2 (ci) 9 sinh {c\z) + ... + 2 (c p ) 9 sinh (c p z) , q > (9) 

C 9 (2) = cosh 2 + 2 (ci) 9 cosh (ci^) + ... + 2 (c p ) q cosh (c p z) , q>0 (10) 
We have: 

^ = (*) ; = (*) (11) 



Using recurrently the formula ({4] 8.486.4): 



and the notation: 



d TL 

—I n {z)--I n {z) = I n+l {z) (12) 



r»W=(~] 0>(*), (13) 



we get: 



-^T n = 4(,)+2,^--J £/ 4pfc (,). (14) 

For n = 1,2,3,4 : 

7\ (z) = z" 1 ^ (2) , T 2 = -z" 3 ^ (2) + z~ 2 C 2 (z) , (15) 

T 3 (z) = Sz' 5 ^ (z) - 3z- A C 2 (z) + z" 3 S 3 (z) , (16) 
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T 4 (z) = -15z- 7 ,Si (z) + 15z- e C 2 (z) - 6z- 5 S 3 (z) + z~ A C A (z) 
The general expressions are: 
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We get: 
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and the following recurrence relations for the coefficients otq : 
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Other general expressions of the coefficients are: 
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0! *" ' 1! " (ra-4)! 
Ignoring the series in the r.h.s. of (Tl4|) . we get approximate expression for I n . Let us give 
here these expressions for the value p = 2 and n — 0, 1, 2, 3 : 



J M (2) = 1.(1 + cosh + 2 cosh , . 

4 V v^/ 
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jM (z) = i. ^sinh^ + v^sinh -^=j , 



- • f — ( sinhz + a/2 sinh —=} + cosh 2 + cosh — = ) , 
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cosh z + cosh — = 1 
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- - ( cosh z + cosh — ^= 1 + | sinh z H -= sinh —= 
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According to the Table 1, even for this very small value of p, the approximation provided 
by these functions, for "moderate" values of the argument (z < 4), is very good. 

Table 1 



z 


l 


2 


3 


4 


I { ap) -Io 
Io 


1.6 x 10" 7 


2.4 x 10" 5 


3.3 x 10" 4 


1.7 x 10" 3 


h 


2.3 x KT 6 


1.4 x 10 -4 


1.2 x 10~ 3 


4.4 x 10^ 3 


I { 2 aP> -I 2 
h 


7.1 x 1(T 5 


10~ 3 


4.5 x 1(T 3 


1.2 x 1(T 2 


h 


1.8 x 1(T 3 


7.3 x IO" 3 


1.7 x 1(T 2 


3 x IO" 2 



It is visible that the precision of the approximation decreases with the order n of the 
Bessel function I n . We can increase it arbitrarly, by increasing the value of p. 

The extension of these results to Bessel functions of real argument is trivial, using the 
formula (El) and: 



S q (-iz) 



[sin z + 2c\ sin (c\z) + ... + 2c q p sin (c p z)~\ , 



(28) 



Cg (—iz) = cos z + 2c\ cos (c\z) 



+ 2c q p cos (c p z) 



(29) 



In conclusion, we have proposed a controlled analytic approximation for Bessel functions 
of integer order. The first 4p — n terms of the series representation of I n is generated exactly 
by the first 4p — n terms of the elementary functions in the l.h.s. of eq.(14). So, our formulae 
can be used, for instance, to find the series expansions of the powers of Bessel functions, a 
subject discussed recently by Bender et al {5] 

The results presented in this paper can be applied to a large variety of problems, mainly 
with cylindrical symmetry, involving Bessel functions at "moderate" arguments. They may 
provide also a useful "visualization" of J n and I n in terms of elementary functions. The 
method cannot be used for asymptotic problems. 
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